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Abstract 

Using dynamical systems methods, we describe the evolution of a 
minimally coupled scalar field and a Friedmann-Lemaitre-Robertson- 
f Walker universe in the context of general relativity, which is relevant 

1} ■ for inflation and late-time quintessence eras. Focussing on the spatially 

flat case, we examine the geometrical structure of the phase space, 
locate the equilibrium points of the system (de Sitter spaces with a 
constant scalar field), study their stability through both a third-order 
perturbation analysis and Lyapunov functions, and discuss the late- 
en . time asymptotics. As we do not specify the scalar field's origin or its 

potential, the results are independent of the high-energy model. 

U ' 

PACS numbers: 98.80.H, 98.80. Cq. 

Keywords: scalar field cosmology, inflation, quintessence, dy- 
namical systems. 

1 Introduction 

The standard cosmological model based on the spatially homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) metric is very suc- 
cessful at describing many observations at different scales. It has become 
normal to include an inflationary epoch [1] in this model during the early uni- 
verse. Although there is no direct proof that inflation actually occurred, and 
other scenarios should still be considered, the 1992 discovery of temperature 
fluctuations in the cosmic microwave background by the COBE satellite [2] 
provided evidence of a nearly scale-invariant spectrum of primordial density 
perturbations of the kind predicted by inflationary scenarios. In addition, 
the study of these temperature fluctuations initiated by COBE ushered in an 
era of "precision cosmology" continued with later cosmic microwave back- 
ground experiments, most notably WMAP and PL A NCK J3l H]- Most mod- 
els of early universe inflation are based on scalar fields, and those based on 



*Physics Department and STAR Research Cluster, Bishop's University, 2600 College 
Street, Sherbrooke, Quebec, Canada JIM 1Z7 

^Physics Department, Bishop's University, 2600 College Street, Sherbrooke, Quebec, 
Canada JIM 1Z7 



1 



quadratic quantum corrections to the Einstein-Hilbert action ( "Starobinsky 
inflation" [5]) can be reduced to the study of scalar field degrees of freedom 

A second revolution in cosmology occurred in 1998 with the discovery, 
obtained by studying type la supernovae, that the current expansion of the 
universe is accelerated [7]. In the context of general relativity, on which 
the standard A-cold dark matter model is based, this acceleration can only 
be explained with a cosmological constant A of extremely fine-tuned, but 
non-vanishing, magnitude, or with a very exotic fluid having pressure P and 
density p related by the equation of state P — p, and dubbed "dark en- 
ergy" . Most models of dark energy are based on a scalar field (ft (also known 
as "quintessence") rolling in a flat section of its potential V(4>). Alternative 
scenarios, seeking to replace the Einstein-Hilbert action ("/(i?)" or "mod- 
ified" gravity [6]), can again be reduced to the dynamics of a scalar field 
degree of freedom. 

Both inflation and quintessence models mandate a general understand- 
ing of scalar field dynamics in general-relativistic cosmology. Furthermore, 
a scalar field provides the simplest field theory of matter, and although 
no fundamental classical scalar field has been discovered in nature so far 
(except possibly for quintessence), they do provide a toy model useful for 
understanding many basic theoretical features of more realistic field theories, 
without the extra details and complications. As such, scalar field theory also 
constitutes an excellent pedagogical tool used in most relativity textbooks. 

In this paper we approach the spatially homogeneous and isotropic cos- 
mology of scalar fields minimally coupled to gravity from the phase space 
point of view. Although dynamical system methods have been widely used 
in cosmology since the 1960s [8] and this type of analysis has been per- 
formed for non-minimally coupled scalar fields [9] and general scalar-tensor 
or f(R) gravity |1U| 111] , we could not find in the literature a complete and 
self-contained analysis for the simpler case of relativity with a minimally 
coupled scalar field, apart from specific scenarios corresponding to particu- 
lar choices of the scalar field potential V(<p), which abound in the literature 
(e.g., [H]-[T5], see also [16] — the literature on specific scenarios is very large 
and here we limit ourselves to quote the papers whose approach is closest 
to the general one that we a dopt)Q By contrast, here we do not commit 
to any particular scenario, and at most, we make general assumptions on 
properties of the potential (such as boundedness or monotonicity), refrain- 
ing from choosing specific forms of the function V ((/)). Given that there is 
no preferred scenario of inflation or quintessence, general considerations are 
valuable. 

1 Ref. [17] studies the minimally coupled case of interest here and presents some of the 
features of scalar field cosmology derived in the following, but its main interest is in initial 
singularities and particular classes of potentials. 
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Since the relevant equations, which reduce to ordinary differential equa- 
tions (ODEs) in this case, are still non-linear and not amenable to exact 
solution, the phase space view becomes important in gaining a qualitative 
understanding of the solutions without actually solving the field equations. 
It is generally believed that in order to say anything about the phase space 
and the qualitative behaviour of the solutions of the equations, one must first 
fully specify the scenario of inflation or quintessence being studied. While 
this is certainly true if one wants a complete qualitative picture of the dy- 
namics, many aspects of the phase space portrait are common to most, if 
not all, scenarios and the study of these aspects, without committing to any 
particular scenario or potential V((j>), is a necessary preliminary for more 
detailed analyses of specific models. The purpose of this paper is to discuss 
these general features, specifically the geometry of the phase space, the exis- 
tence, nature, and stability of the fixed points, and the late-time behaviour 
of the solutions, without specifying the form of the scalar field potential en- 
ergy density, and instead making some generic assumptions on its behaviour 
(boundedness, presence of asymptotes, etc.). 



2 Background 

We consider a scalar field minimally coupled to the spacetime curvature as 
the only source of gravity in the Einstein field equations. This assumption 
is fully justified in inflationary scenarios of the early universe |15| . and only 
approximately justified in quintessence models of the late universe |18j . In 
the latter case, scalar field dark energy is present along with a dust fluid, 
which combine to determine the dynamics of the universe. However, ob- 
servations suggest that dark energy comes to dominate the dynamics very 
quickly, starting from redshifts z ~ 0.5, thus we can once again neglect the 
dust fluid and other forms of energy in the late regimes. In short, there is 
plenty of motivation to study scalar field cosmology. 

The Lagrangian density of a scalar field <j) minimally coupled to the 
spacetime curvature is^| 

£ W = _l V a 0V a <A-y(0), (2.1) 

where V{4>) is the scalar field potential. The action for gravity and the scalar 
field is 



S = <Tx 



+ £ W (0, 9^) 



16vrG 



S(g) + <%) , (2.2) 



where g^ v is the spacetime metric, g is its determinant, and R is its Ricci 
scalar. The action (|2.2p is also the action for general scalar-tensor gravity 



2 We follow the notations of Ref. [H]. 
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in vacuo, after performing a conformal transformation to the Einstein frame 
(e.ff., [20]). The variation of the scalar field action Su) = f ofx^J^g 
gives the stress-energy tensor 

_9 1 
T$ = ^|^7 = V M 0V^ - - ^V>V a </> - V{<j>) . (2.3) 

A spatially homogeneous and isotropic universe is described by the FLRW 
line element 



ds 2 = -alt 2 + a 2 (t) 



-^—+r 2 (^ 2 + sin 2 0V) 
1 — kr A 



(2.4) 



in comoving coordinates (t,r,0,ip), where a(t) is the scale factor and k is 
the curvature index. The Einstein equations 

R^u-\g^R = ^GT IMU (2.5) 

(where R^ u is the Ricci tensor and R = -R%) reduce to ODEs for the scale 
factor and matter degrees of freedom. It is customary to approximate the 
matter content of the universe with a single perfect fluid with four-velocity 
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§0/i j n C omoving coordinates, energy density p, pressure P, and energy- 
momentum tensor 

Tfo, = (P + p) u^u v + P 9fMV . (2.6) 

The pressure and energy density are usually related by a barotropic equation 
of state P = P(p), often of the form P = wp where the constant w is called 
the "equation of state parameter". The Einstein equations (|2.5|) in the 
presence of a single perfect fluid reduce to 

+ (2.7) 

a d 

a 2 8ttG k . . 

P--5, (2-8) 



a 2 n2 



p + 3-(P + p) = 0, (2.9) 



a 



where an overdot denotes differentiation with respect to the comoving time t. 
Eqs. (|2.7|) and (|2.8p are called the acceleration equation and the Hamiltonian 
constraint, respectively, and the Klein-Gordon equation (j2.9|) is nothing but 
the covariant conservation equation V V T^ U = (when <f> ^const.). The 
Klein-Gordon equation is not independent of eqs. (|2.7p and fj2.8[) and can 
be derived from them. Excellent pedagogical analyses of the phase space of 
a FLRW universe coupled to a perfect fluid are available in the literature 

eh. 
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In a FLRW universe, a gravitating scalar field must necessarily depend 
only on the comoving time, <j> = <f>(t), in order to respect the spacetime 
symmetries. Therefore, its gradient is timelike (or null but trivial if 
4> =const.). In regions where X7 a 4>V a (ft < 0, we can introduce the four-vector 

u = V ^ (2.10) 

Viv^v^i 

with u^u^ = —1, and the scalar field is equivalent to a perfect fluid with 
stress-energy tensor of the form (|2.6|) and energy density and pressure |22] 

p = y + v W> (2 - n) 

P = ^--V(4>). (2.12) 
One can define the effective equation of state parameter 



MM)s i^pjm (2 . 13) 



The Einstein-Friedmann equations f|2.T[) — fj2.9[) become 

a _ 8vrG 
a 



<]> 2 - V {<!>)) , (2.14) 



a\ 2 8itG 



j +V(</>)\ --, (2.15) 




+ 3-0+ ^- = 0. (2.16) 



a 

In the following it will be useful to rewrite these equations in terms of the 
Hubble parameter H = a/a as 

H = -H 2 -^{tf-V^))+- = -ATtGtf + - 2 , (2.17) 

a z a z 




H 2 = + (2.18) 



4> + 3H(f> + V =0, (2.19) 

where a prime denotes dfferentiation with respect to <f>. These equations can 
also be derived from an effective Lagrangian or Hamiltonian [23] , 

The equations of scalar field cosmology are non-linear and few exact 
solutions are known for particular choices of the potential V{<j>). We would 
like to discuss the dynamics of the variables a(t) and 4>(t) in as much depth 
as possible without choosing a specific form of V{4>). Before we begin, let 
us note that 
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• for V((f>) = the scalar field is equivalent to a fluid with stiff equation 
of state P = p, which does not seem to be very relevant for inflation 
and late-time acceleration (although it is relevant for matter at nuclear 
densities in the core of neutron stars, and possibly near the Big Bang 
singularity [17]). 

• For V{4>) = Vo = const, the potential reduces to a pure cosmological 
constant A. The scalar field stress-energy tensor (|2.3p reduces to 

T^ v = ~^q9^u ~ d^4>d v (j) - i g^ u d a (f)d a <j) , (2.20) 

with A = 8ttGVo. Further, for <j> = (j)Q =const., one recovers the stress- 
energy tensor of a pure cosmological constant. 

3 Phase space 

Eqs. (|2.14p and (|2.16p describe the evolution of a(t) and 4>(t) (remember 
that there are only two independent equations in the set (|2.14|) - (|2.16p if 
4> is not constant). Eq. ()2.15p is a first order constraint (contrary to eqs. 
(|2.14p and (|2.16p which are of second order). The phase space is, therefore, 

a 4-dimensional space ^a, d, (ft, <fj , but the Hamiltonian constraint (|2.15p 
forces the orbits of the solutions to live on a 3-dimensional hypersurface, 
introducing a relation between the four variables. For example, one can 
use the constraint to express <fi in terms of the other three variables, 4> = 
<j>(a,a,</>). 

For particular choices of the scalar field potential, and especially for 
k 0, one can change variables to functions of a, a, (ft, (/) which can lead to 
exact solutions or to simpler calculations. In general, however, these new 
variables do not have an immediate or clear physical meaning and are to be 
regarded as a mere mathematical trick to perform calculations. Often the 
results of these calculations cannot be translated explicitly or easily in terms 
of the variables ^a, d, <fi, <pj . However, current observations seem to indicate 

that we live in a spatially flat (k = 0) universe, which is much simpler to 
analyze than the k 7^ case. This is the situation that we consider in the 
following. 

4 Spatially flat FLRW scalar field cosmologies 

The description of the phase space greatly simplifies for k = as, in this case, 
the scale factor a(t) appears in the dynamical equations only through the 
combination a/a = H, the Hubble parameter, which is a physical observable 
obtained by fitting theoretical models to cosmological data. Since <fi is the 
only matter field in the theory, it is natural from the field theory point of 
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view to choose it as another dynamical variableO By choosing H and <p 
as dynamical variables, the phase space reduces to the 3-dimensional space 
[H, (f>, <£) , but the orbits of the solutions of eqs. <^JJ!)-^Wj with k = 
are forced to move on a 2-dimensional subset of the phase space by the 
Hamiltonian constraint (|2.18p . 

Let us examine the structure of the "energy surface" on which the orbits 
are forced to move. We choose to eliminate (f> by expressing it in terms of the 
other variables (H, <p) in eq. (|2.18p with k = 0, which can then be viewed 
formally as a quadratic algebraic equation for cp and solved, obtaining 

For certain choices of the potential V{<p), an arbitrary choice of values of the 
pair (H, (p) could make the argument of the square root on the right hand 
side negative. Therefore, in general, there can be a region of the phase space 
forbidden to the orbits of the dynamical solutions, 

JF= {(H, : 3H 2 < 8irGV((p)} (4.2) 

("forbidden region"). This region may or may not exist depending on the 
form of V (cj)). 

There are two portions of the phase space region accessible to the dynam- 
ics (the "energy surface" corresponding to vanishing effective Hamiltonian 
[23]), corresponding to the two signs of the right hand side of eq. (|4.1D . 
These sets are symmetric with respect to the <p = plane of the (^H, cp, cpj 
space. We call these two subsets of the energy surface "upper sheet" and 
"lower sheet" , corresponding to the positive and negative sign, respectively. 
In the upper sheet (p is always increasing i^cp > 0^ , while on the lower sheet 

<p is always decreasing \ 6 < 0^. The two sheets are either disconnected, or 
always join on the plane <p = 0, on the boundary of the forbidden region 

B = dF = { (H, (P, $\ : <p = ^ 3H 2 = 8ttGV((P)} . (4.3) 

Figs. [H an d [5] show the upper and lower sheet for the example potential 
V{cp) = m 2 (j) 2 /2. The dynamics of the spatially curved (k ^ 0) scalar field 
universe are confined to either side of the "energy surface" corresponding 
to k = in the phase space — this fact can be deduced by reducing the 
Hamiltonian constraint to 

3 An astronomer would instead choose the density of the matter field fi^ (in units of 
the critical density) as another variable. 



7 



Figure 1: The upper sheet corresponding to the positive sign in eq. (|4.ip . 
for the quadratic potential V((f>) = m 2 (p 2 /2 (in arbitrary units). 

Trajectories corresponding to k > would exist above the k = upper sheet 
(i.e., for larger values of (p than those corresponding to the k = upper 
sheet) and below the k = lower sheet (i.e., for lower values of <ft). Trajec- 
tories corresponding to k < would exist between each k = sheet (i.e., for 
values of <j) comprised between those given by eq. (|4.ip ). This property was 
realized in Ref. [p3] for the specific inflationary potential V = m' 2 (j) 2 /2, but 
the conclusion is general. Trajectories in a region corresponding to k > 
cannot cross the k = sheet and move to regions corresponding to k < 0, 
and vice-versa. Such dynamical transitions betwen different topologies of 
the universe are forbidden (the topology of spacetime is not ruled by the 
dynamics). 

Finally, the lower dimension of the "energy surface" leads one to believe 
that chaos is impossible in the dynamical system under study. This state- 
ment is not trivial given that the standard results on the absence of chaos 
in a two-dimensional phase space are proven for a plane, not for a curved 
surface or for a subset of a higher-dimensional phase space obtained by glu- 
ing two 2-dimensional sheets [M]- However, it is not difficult to reduce this 
situation to the standard C3jSG ? clS has been shown for scalar-tensor gravity 
in [2S] . The theory of a minimally coupled scalar field in Einstein gravity is 
contained in this reference as a special case. 
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Figure 2: The surface describing the Hamiltonian constraint (|4.ip for the 
quadratic potential V(cj)) = m 2 (j) 2 /2 (in arbitrary units). The upper and 
lower sheets join at the <p = plane to form a cone. 

5 Equilibrium points 

Having chosen H and <f> as dynamical variables, the equilibrium points of the 
dynamical system (when they exist) are, by definition, of the form (^H, tfij = 

(0,0) and yH,(pj = (0,0), or {H,<p) = (Ho,4>q) = (const., const.), and they 

must all lie in the eft = plane, and therefore, on the boundary B of the 
forbidden region (if this region exists). These equilibrium points are de 
Sitter spaces with a constant scalar field. When they exist, they are the 
only de Sitter spaces possible in this theory. In fact, eq. (|2.17|) with k = 
reduces to H = —AttG4> 2 , and a de Sitter space with H =const. necessarily 
has 4> = const, as wellQ A degenerate case is Hq = 0, which corresponds 
to Minkowski space, de Sitter spaces are important in cosmology because 
they are usually attractors in inflation and quintessence models [151 EE]- 
For <j> =const., and reduce to = -V = -V{cj) ) and T$ = 
—Vq g^ u , i.e., to a pure cosmological constant term with A = SttGVq. 

The necessary and sufficient conditions for the existence of de Sitter fixed 



4 By contrast, with non-minimally coupled scalar fields, de Sitter spaces with a non- 
constant scalar field are possible [S]. 
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Figure 3: Trajectories converging to a Minkowski fixed point \^H, (ft, (ft 
(0,0,0) for the previous example V{4>) = m 2 (ft 2 /2. 
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points are easily obtained from eqs. (j2. 17[) (j2. 19[) with k = 0: 
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V , 



(5.1) 



3 








(5.2) 



which obviously require Vq > (Minkowski space is obtained for Vo = 0). 
Eq. (|5.2p expresses the condition that V{(ft) has an extremum or a point 
with horizontal tangent at </>o- 

Fig. [3] shows two trajectories, corresponding to different initial condi- 
tions, converging to a Minkowski space attractor point for the example of 
the V((ft) = m 2 (ft 2 /2 potential. 

Attractors (or repellors) could exist as an asymptotic limit in a or (ft. To 
check for these we must search for fixed points with infinite values of the 
variables. 

5.1 Fixed points at infinity with a Poincare projection 

Fixed points at infinity can be found by adopting polar coordinates (r, 9) 
with 



H = r cos 9 



= r sin 



9 



(5.3) 



and the standard Poincare transformation r — > r with 



r 




I 



(5.4) 
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which maps infinity onto the circle of radius r = 1. Since 



H 



(1 + r) j_ 

r cos I 



r (1 — rY 
1 + r) 



r sin # + 



1 — r 



1-r 



' sm( 



' COS I 



(5.5) 



(5.6) 



2Vr (1 - rf 

fixed points (H,4>) = (const., const.) correspond to (r,8) = (const., const.) 
thanks to the linear independence of the sine and cosine functions. The 
dynamical system (|2.17p - (l2.19p becomes 



(1-r) 2 



4vrG 



cos 2 



/ r ^ + r) 
\2y/¥(l-r] 



sin# + 



1 — r 



cos 0) + 2V 



, (5.7) 



r (1 + r) 

— ^ =■ COS (/ — ■ 

rfl-r) 

F(l + r) 



1 sm( 



-4vrG 



-, 2 



2y¥(l 



-\2 



sm( 



1-r 



' COS I 



(5.8) 



? 

j iy*ty* j 



2yr (1 



_l2 , ^2_ j_2 , i_2_ 

ry I iyr ryt ty _l_ y« 

4 F 3/2 ^ _ f )2 + 4 ^ ^ _ -y 



sm( 



+ —je zn e cos 6 + 

yr (1 — r) 
r 



1 — r 



1 cos l 



sm( 



1 cos" + 



r (1 + r) 

.(l-r) z 2 (1-r) 3 

Setting (H,(j)) = (0, 0) and (h, <j>) = (0, 0) yields 



sin I 



V 



r cos 



(1-r)' 



Vr 







0. 



(5.9) 



(5.10) 



(5.11) 



where 4>q = (f>(ro, 8o)- In order to satisfy eq. (|5.1U|) in the limit r — > 1 we 
must have either 

1. cos6> = 0, corresponding to H — > 0, — > dboo, and V — > dboo) = 
(this situation includes potentials V(<f>) with compact support). 
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2. cos# = ±1, corresponding to H — > ±00, — > 0, and V (0 — > 0) =00 
(i.e., V has a vertical asymptote at = 0. This is the case of the 
potentials V(0) oc 1/0°, a > used in many quintessence models). 

3. 9 7^ 0, ±7r, ±7r/2, which allows -ff — >■ ±00, — > ±00, and V (0 — )■ ±00) = 
00. This case includes unbounded monotonic potentials such as V((p) = 
Vo e ±a< ^ (scalar field cosmology with exponential potentials is studied 
in detail in Ref. [H]). 

Fixed points corresponding to any of these situations must have a po- 
tential that asymptotically satisfies eq. (|5.1ip as well as the stated condi- 
tions. A fixed point satisfying the conditions of situation 1) corresponds to 
Minkowski space with no potential. Situations 2) and 3) both correspond to 
extreme cases of de Sitter space, with their potentials diverging. A possible 
situation corresponding to case 3) is V((p) = Vo In (0/0o). 

The next question that one can ask regards the stability of these equi- 
librium points. 

5.2 Stability with respect to homogeneous perturbations 

It seems intuitive that if V(4>) has a local minimum at 0o, and a de Sitter 
equilibrium point (Hq^q) satisfying eqs. (|5.ip and (|5.2p exists, it will be 
stable, and vice-versa, it will be unstable if V((j>) has a local maximum. 
However, this statement would be a bit naive because <p(t) couples to the 
other variable H(t) and one must consider both variables simultaneously. 
Here we consider homogeneous perturbations of the fixed point (Hq,4>q), 
i.e., we write 

H(t) = H + e 8 x H(t) + e 2 <5 2 ^(t) + e 3 5 3 H(t) + ... , (5.12) 

0(t) = 0o + e 5 1( />{t) + e 2 5 2 m + eH^(t) + ... , (5.13) 

where e is a smallness parameter and 5^H and 5^4> depend only on time. 
In general, one should consider more general perturbations 5H(t,x), S(f>(t,x) 
and even anisotropic perturbations. Inhomogeneous perturbations are sub- 
ject to notorious gauge dependence problems and can only be treated rig- 
orously in the context of a gauge-invariant formalism [26]. This kind of 
formalism is necessarily very detailed and complicated and the correspond- 
ing gauge-invariant variables are susceptible to physical interpretation only 
after a gauge is fixed (and then, different gauges produce different interpre- 
tations). For clarity, and to avoid the high degree of sophistication needed, 
we will confine our analysis to homogeneous perturbations. Therefore, if a 
de Sitter fixed point is stable with respect to homogeneous perturbations it 
may still be unstable with respect to inhomogeneous ones. In the following 
we assume that an equilibrium point exists. 
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By inserting eqs. (pTT2|) and (fSTTB"]) into eqs. (f2T71) - ([2J9|) with k = and 
using the zero order eqs. (|5.ip and (|5.2p for the equilibrium point (Ho,(f>o), 
one obtains the perturbed Hamiltonian constraint 



2eH 5 1 H + 2e z H 5 2 H + e 2 5 1 H 2 + 2e 4 Ho5 3 H + 2e 6 5iH5 2 H 



8ttG 



(5.14) 



the acceleration equation 

e<5iij + e 2 J 2 ii~ + e 3 5 3 H = - (2eH 8 1 H + 2e 2 H 5 2 H + 2e 3 H 6 3 H 



+2e 3 5 1 H5 2 H + e 2 5 1 H 2 J 

and the Klein- Gordon equation 
€5x4) + e 2 5 2 (f> + e 3 5 3 (f> 



8irG 



e 2 5 1 $ 2 + 2e 3 5 1 & 2 c, 



Vo (e 3 5 1( j> 3 ) 



(5.15) 



-3 (eHoSnj) + e 2 H 5 2 4> + e 2 5 l H5 l 4> + e 3 H 5 3 <j) + e 3 5 1 H5 2 ^ + e 3 5 2 H5 l ^ 



In 



+V (e<5i0 + e 2 <5 2 + e 3 S 3 cf) + ( e 2 <5 1( /> 2 + 2e 3 5^ 5 2 



Vr 



(IV) 



Jr J.3 



6 



0, 



(5.16) 



where V = V (cpo)i e ^ c - ^(<A) an d dV/d(f) have been expanded to third 
order as 

V" V'" 
^ (0) = Vo + -f [z&i<Kt) + e 2 5 2 <P(t) + . . . ] 2 + ^- [edx^t) + . . . f + . . . , 

(5.17) 



V\<j>) = Vo [e5 1( f>(t) + e 2 5 2 <P(t) + e 3 5 3 <j ) (t) + . . . ] 



V'" o 1/ (IV) 

[ e ^(t) + e 2 W)+ ...] 2 + ^[e5 1( A(t)+ ...] 3 + ... , 



(5.18) 
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using the fact that Vq = 0. To first order in e, eq. (|5,14p yields 2Hq5\H = 
and, if the de Sitter equilibrium point is not a degenerate Minkowski space 
with H o = 0, therU 

5iH = (5.19) 



and eq. (|5.15p then yields 8\H = 0. The Klein-Gordon equation (|5.16|) then 
decouples and reduces to the equation of the damped harmonic oscillator 



5 1 4>+m Q 5 1 (i) + y o <5i0 = o. 



(5.20) 



a second order ODE with constant coefficients. The associated algebraic 
equation is 



A + 3F A + V 



0. 



which has the roots 



Ai, 



3Hn 



AVI 

ml 



(5.21) 



(5.22) 



If 9-ffg - 4Vq / 0, the general solution of eq. (fOUj) is 



SMt) 



-3H t/2 



+ C 2 exp 



Ci exp 



<9H*-W " 



9flg - 4y " 



2 



Cie Alt + C 2 e A2i , (5.23) 



where Ci 2 are arbitrary integration constants. It is easy to see that V > 
is required for stability. In fact, if Vq > 0, then 9Hq - AV " < 9Hq. 

If V " > 9^/4 = 67rGV , then X lj2 = 



-3H ±i^\9H$-4V "\) /2 
have imaginary parts and the two independent modes are of the form 



exp 



-SHpt \ 



■ exp 



±ij 9F 2 -4F " 



(5.24) 



which decay because of the first exponential if Hq > 0, and increase without 
bound if Hq < 0. 

If Vq = 9iZg/4 then A x = A 2 = -3H /2 and the solutions of eq. ([5T20D 

are 



(t) 



-3H t/2 



{Ci + C 2 t) 



(5.25) 



It is an old adage in cosmological perturbation theory that there are no linear pertur- 
bations of de Sitter space sourced by a (minimally coupled) scalar field [27] . 
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again, stable if Hq > 0. 

If < V " < 9^/4, then < 9H% -AV " < 9H$, or ^9i^ - 4y o " < 3iJ 

and Ai 5 2 < 0, so that the independent modes e Al - 2 * do not grow. 
Finally, ifV =0 in eq. ()5.20p . the solution is 

SKf>(t) = C 1 + C 2 e- 3Hot ; (5.26) 

this solution is stable for Hq > 0, but not asymptotically stable, as it does 
not fully decay. We can therefore conclude that: 

• If Hq > 0, then the de Sitter equilibrium point (Hq,4>q) is asymp- 
totically stable if V > 0, unstable if V < 0, and stable but not 
asymptotically stable if Vo = 0. 

• If Hq < 0, the de Sitter equilibrium point (Hq,(/)q) is always unstable 
(to first order and all higher orders). 

Furthermore, to first order there is no perturbation 5\H and the perturba- 
tions 5H and 5(f) are decoupled. 

The second order equations yield 

hR = If + y °''^ 2 ) ' (5 ' 27) 

5 2 H = -4vrG5i0 2 , (5.28) 

V" 

5 2 4> + 3H 5 2 <j) + Vq5 2 $ = — ^-5i0 2 . (5.29) 

The second order perturbations 5 2 H and 6 2 (f> depend on the first order ones 
5\4> and their derivatives. These act as a source term in the Klein-Gordon 
equation (|5.29|) for 5 2 (f), which is a damped forced harmonic oscillator equa- 
tion. From the general theory of ODEs, it is clear that if Hq > 0, the friction 
term 3Hq5 2 4> will correspond to positive friction, while if Hq < 0, there is 
"anti-friction" leading to instability. Therefore, 

• For Hq > 0, there is stability if Vq > 0, and instability if < 0. 

• If Hq < 0, the perturbations 5 2 (p grow without bound and the de Sitter 
equilibrium point (Hq,4>q) is unstable. 
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To third order in e, one obtains 
4ttG 



5*H 



6 1 <j>5 2 <t> + VQ5 1 <i>5i<i>+ ^Oi^ 3 



-87rG$i0$ 2 0, 



(5.30) 



(5.31) 



6 3 (t> + 3H Q 6 3 <I> + V Q 53<, 



2ttG 
3^ 



Vr 



(IV) 



6 



$i< 



(5.32) 



By considering only situations in which we have stability at the lower orders, 
the lower order perturbations and their derivatives, which are acting as 
sources in the higher order equations, are bounded, and eqs. (|5.30p and 
(|5.3ip guarantee that 5 3 H and its derivatives are bounded by (small) initial 
conditions. Eq. (|5.32p again takes the form of a driven, damped harmonic 
oscillator with instability if Hq < 0, or if Hq > with V < 0, and stability 
for H > and Vq > 0. 

There remains the case of the Minkowski fixed point (Hq, 0q) = (0, 4>q), 
a degenerate de Sitter space. This situation occurs if Vb = and V Q = 0. 
Eqs. (|27T7|) - pn$|) then reduce, to third order, to 

e 2 5 1 H 2 + 2e 3 8 1 HS 2 H 



8ttG f 1 
"3" I 2 



e 2 rf + 2e 3 $!# 2 0) 



+ \ Vq {e 2 ^ 2 + 2e 3 5^5 2 4>) + ± V "' (e 3 5 1( f> 3 ) \ , 
edxH + e 2 6 2 H + e 3 5 3 H = -2e 3 5iH5 2 H - e 2 6 1 H 2 



(5.33) 



8ttG 



e 2 8 1 <j) 2 + 2e 3 5!<j>5 2 c 



V (e% 



3x J.3> 



(5.34) 
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(e6i4> + e 2 5 2 '4> + e 3 5 3 ^j 

+3 (e 2 5 1 H5 1 <p + e 3 5 1 H5 2 <p + e 3 5 2 H5 1 c 

,,, 

+Vq (e<5i0 + e 2 b 2 <\> + e 3 5 3 (f) + -2- (e 2 5^ 2 + 2e 3 '5^ 5 2 <f) 

+^ — e 3 5 1 cp 3 = 0. (5.35) 
6 

The friction term is now absent in eq. (|5.35|) . To first order one obtains 
8\H = and 

5i4> + Vq 'Sxt = , (5.36) 

therefore there is stability for Vq > 0. 
To second order we have 

^^[ht^y^. (5 . 37) 

(which yields no information on second order perturbations), 

S 2 H + dxH 2 = ^ (-S^ 2 + i F "<^ 2 ) = -inGd^ 2 , (5.38) 

V" 

5 2 4> + Vq5 2 (I) = -36xH6ij> ■ ( 5 - 39 ) 

Eq. (|5.38|) guarantees stability in 5 2 H, while eq. (|5.39p guarantees stability 
to second order if V > 0. 
To third order we have 

S 2 H = (siM + V^5 1 <j ) 5 2 (t> + 1 F "'M 3 J , (5.40) 

assuming <5iff 7^ (again, this equation provides no information on third 
order perturbations), 

5 3 H = -28 X H8 2 H - ^ Us x M - F o %# 2 - i Oi</> 3 ) , (5-41) 



/ ta('V) \ 

5 3 + <<5 3 = - 35^5^ + 3£ 2 #<M + Vo" 6 1( j)5 2 (t) + <5i0 3 , 

(5.42) 

which again provides stability for V > 0. 
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5.3 Inhomogeneous perturbations 



A treatment with respect to inhomogeneous perturbations is necessarily 
more complicated, requiring the use of a gauge-invariant formalism (see, 
e.g., |26j). To first order (for which the gauge-invariant formalisms apply), 
the results on the stability of de Sitter spaces already obtained also hold for 
inhomogeneous perturbations. In fact, the stability of de Sitter spaces in 
the very general theory of gravity described by the action 



S 



f{R,4>) u{4>) 



g^V^Vrf - V(<f>) 



(5.43) 



was studied in [25]. This class of theories contains scalar-tensor gravity if 
f(R, 4>) = f(4>)R and higher order gravity if f(R, <fr) = f(R) and <fr =const. 
A de Sitter space (iJo,<^o) is stable with respect to inhomogeneous pertur- 
bations if and only if [28] 



11- 



fa 



1 + 



< 0. 



(5.44) 



General relativity with a minimally coupled scalar field corresponds to the 
trivial case 

/(0,jR) = 8^' W = (5 ' 45) 

which yields the first order stability condition V " > 0. This result is gauge- 
invariant and reproduces, to first order, the one which we have already 
obtained for homogeneous perturbations (the equivalence between homoge- 
neous and inhomogeneous perturbations in this respect does not extend to 
FLRW universes other than de Sitter space). 



5.4 Lyapunov functions 

The Lyapunov method [24] can be easily applied to scalar field cosmology 
in order to assess stability non-perturbatively, and to estimate the size of 
the attraction basins of stable fixed points in the phase space. 

Let (Hq, 4>q) be a fixed point of the dynamical system fj2.17|) f)2.19[) : then, 
if V(cfr) has a local minimum at cj)Q, and Hq > 0, the C 1 function 

Li (h,</>,£) = y+^-Vb (5.46) 
is a Lyapunov function. In fact, 

• L\ (H, 4>,A>0 in a domain T> containing (fro, except at (Hq, (fro); 
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L t (#0,^0,0) = 0; 



• dLi/dt = cf)yj) + V j = —3H(j) 2 is strictly negative in T>, except at the 
fixed point, where dL\jdt vanishes. 

Therefore, the fixed point (iJo,0o) is asymptotically stable. If V(4>) has 
only one minimum at c/)q, the attraction basin of (Hq, 4>q) is the entire phase 
space, while if there are other minima at (f>i,(f>2, ■ ■ • , the attraction basin of 
(Hq^q) is finite and will be limited by separatrices between the attraction 
basins of other fixed points. 

If Hq > and V(i^) has a local maximum at 4>q (with V = 0, Vq < 0), 
then the C 1 function 

L 2 (h, 0, fj = - V(4>) + V (5.47) 

is such that L 2 ^-ff , 0, 0^ < in a domain containing (Hq,(/)q) at which L 2 

vanishes, and dL 2 /dt = 3H(j) 2 > except at the fixed point itself where 
dL 2 /dt vanishes. This guarantees that the fixed point is unstable and is a 
repellor. 



6 Late-time behaviour of the solutions 

Some conclusions on the asymptotic behaviour of the solutions at late times 
can be reached under certain assumptions on the scalar field potential, with- 
out fully specifying the form of V(4>). 

First, for spatially flat universes, eq. (|2.17p implies that H < 0, with the 
equality being satisfied only for the de Sitter fixed points. Hence, outside 
of fixed points, H(t) is always a decreasing function. Assuming that H 
starts out positive and that V{4>) > (which guarantees that the energy 
density is positive), the Hamiltonian constraint f|2. 18[) shows that H cannot 
become negative because, due to continuity, it would have to vanish first and 
the trajectory of the solution in phase space would then cross a fixed point, 
which is impossible. Therefore, if H starts out positive [negative], it remains 
positive [negative]. Let us consider, for definiteness, the case H > 0. Since 
H is bounded from below by zero and H < 0, the graph of H(t) cannot 
cross the H = axis and we must have H > (assuming that H £ C 2 ). 
Since H(t) = —AirGcf) 2 , as H(t) — > at late times so also must (j)(t) — > 0, 
or (H (t) , <j)(t)) — > (_ffo,</>o) (with Hq > 0) at late times. The trajectories 
must asymptotically approach a de Sitter attractor point. One possibility 
is that Hq = 0, corresponding to Minkowski space. This could be a point 
at infinity, in which case the energy content of the universe gets diluted 
in the future expansion and the universe more and more resembles empty 
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Minkowski space. In general, for a strictly monotonic potential which is non- 
negative (or otherwise bounded from below), we always have dV/d<f) 7^ 
and there cannot be equilibrium points. In this case \(p\ — > +00 at late 
times with </> — > due to friction (in an expanding universe), and then 
also H = -AirGcj) 2 -)• 0. The phase space orbit tends to a de Sitter point 
(Ho,(fro) = (87rGVo/3, ±00), where Vq is the asymptotic value of V(4>). 

A further consequence of the acceleration equation H = —4irG(p 2 is that 
there cannot be limit cycles as H and <p would have to repeat themselves for 
periodic orbits, while here H is monotonically decreasing and thus cannot 
be periodic. 

The late-time asymptotics become much more complicated if the scalar 
field couples non-minimally to the Ricci curvature R or if it is a phantom 
field with the "wrong" sign of the kinetic energy [9j [29l [30] . 

7 Conclusions 

It is possible to partially analyze the dynamics of a minimally coupled 
scalar field in FLRW cosmology, described by the coupled Friedmann-Klein- 
Gordon equations (|2. 18() , (|2.17p , and (|2.19p constituting a dynamical system. 
It has been shown in the previous sections how the study of the geometrical 
structure of the phase space and of the stationary points and their attrac- 
tion/repulsion basins can be carried out without specifying the form of the 
scalar field potential V((f>). Limiting the analysis to spatially flat universes 
reduces the dimension of the phase space from four to three. Further, the 
trajectories of the system lie on two intersecting energy "sheets" (corre- 
sponding to setting an effective Hamiltonian equal to zero) given by the 
Hamiltonian constraint (|4.1|) . The two sheets correspond to <j) > and 
(f> < 0. There is no possibility of chaos in this space due to its lower di- 
mensionality. The dynamics of non-spatially flat universes are described by 
orbits constrained to take place above the upper sheet and below the lower 
sheet ^|0| > \4>fiat\^ f° r k > 0, and between these sheets for k < 0. 

Any stationary points are necessarily de Sitter universes with a constant 
scalar field, which includes Minkowski space as the special case Hq = 0. 
These points were found to be stable for Hq > and V > 0, by both third 
order homogeneous perturbation analysis and Lyapunov's second method. 
The size of their attraction basins is dictated by the form of the potential. 
The attraction basins are global if Vq = corresponds to a single sink. The 
possibility of limit cycles is excluded. 

In an asymptotic analysis, if we assume that H starts off positive (an 
initially expanding universe) and V{<j>) > (enforcing the weak energy condi- 
tion), then H (t) will remain positive for all times and all trajectories meeting 
these conditions will be asymptotically attracted to de Sitter spaces. Phys- 
ically, this means that for these conditions the universe will always expand, 
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with the scale factor a(t) becoming exponential at late times. Stationary 
points could exist as an infinite limit in H and/or (f> for certain potentials. 

From the previous results one can see that an FLRW universe contain- 
ing a single scalar field in general relativity has relatively simple dynamics. 
The rather straightforward discussion presented here needs, of course, to be 
supplemented by a more detailed study which can only be carried out by 
fully specifying the potential V(4>), and this is the limitation of the present 
paper. At the same time, not specifying the potential makes our discussion 
completely general. An added bonus of the material presented here is its 
pedagogical value for a general introduction to inflation and quintessence 
models. 
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